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Abstract 

Sustainably managing plant resistance to epidemic pathogens implies controlling 
the genetic and demographic changes in pathogen populations faced with resis- 
tant hosts. Resistance management thus depends upon the dynamics of local 
adaptation, mainly driven by the balance between selection and gene flow. This 
dynamics is best investigated with populations from locally dominant hosts in 
islands with long histories of local selection. We used the unique case of the 
potato late blight pathosystem on Jersey, where a monoculture of potato cultivar 
'Jersey Royal' has been in place for over a century. We also sampled populations 
from the coasts of Brittany and Normandy, as likely sources for gene flow. The 
isolation by distance pattern and the absence of genetic differentiation between 
Jersey and the closest French sites revealed gene flow at that spatial scale. Micro- 
satellite allele frequencies revealed no evidence of recombination in the popula- 
tions, but admixture of two genotypic clusters. No local adaptation in Jersey was 
detected from pathogenicity tests on Jersey Royal and on French cultivars. These 
data suggest that long-distance gene flow (~ 50/100 km) prevents local adapta- 
tion in Jersey despite a century of local selection by a single host cultivar and 
emphasize the need for regional rather than local management of resistance gene 
deployment. 



Introduction 

From an evolutionary ecology point of view, local adapta- 
tion (or maladaptation) between hosts and pathogens 
results from the balance between evolutionary forces acting 
simultaneously on pathogen life history traits involved in 
trade-offs (Gandon and Day 2009). Therefore, understand- 
ing the evolutionary mechanisms acting upon pathogen 
populations, and quantifying their balance, is essential to 
develop durable control strategies of these pathogens. This 
is the basis of Darwinian medicine in humans and animals 
(Nesse et al. 2010; Nesse 2011) and has the same applica- 
tions for crop protection (Andrivon 2012). 

However, predicting local adaptation remains difficult 
for both qualitative and quantitative pathogenicity traits. 
In the first case, corresponding to the boom-and-bust 
cycles observed in gene-for-gene systems, local adaptation 



is basically driven by disruptive selection following the 
appearance of virulent variants in an avirulent {i.e., incom- 
patible) population, through mutation or gene flow (Kiyos- 
awa et al. 1995; Brown and Tellier 2011). The virulent 
mutants are usually assumed to bear a fitness cost over 
avirulent genotypes on the hosts devoid of the matching 
resistance genes, although stable polymorphism can exist in 
the absence of such a cost (Leonard 1994; Tellier and 
Brown 2011). It has been possible to measure fitness costs 
of virulence in a few cases, although these costs are not 
always very high (Leach et al. 2001; Bahri et al. 2009; Mon- 
tarry et al. 2010a), and can sometimes be compensated by 
additional mutations (Lannou 2012). Avirulent genotypes 
will therefore usually persist in the pathogen population 
as long as the selection pressure exerted by the hosts is 
low {i.e., resistant hosts occupy a low fraction of the 
host acreage) or gene flow continuously brings into local 
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populations genotypes maladapted to the local resistant 
hosts. This has been shown to occur in several Blumeria 
graminis populations in Europe (Andrivon and Limpert 
1992). The same evolutionary processes probably apply also 
to quantitative resistance, but the complex nature of this 
type of resistance makes a prediction all the more difficult. 
Indeed, quantitative resistance often results from the alter- 
ation by the host of several pathogen life history traits, such 
as latent period duration, host colonization rate, or repro- 
ductive capacity (Clement et al. 2010). It is therefore 
polyphenic and polygenic, which makes an assessment of 
its fitness impact on pathogen (and host) genotypes both 
theoretically and experimentally challenging (Pariaud et al. 
2009; Montarry et al. 2010b; Lannou 2012). 

Theoretical models can be used to predict the occurrence 
of local adaptation (or lack thereof) in pathogen popula- 
tions. They are based either on the comparisons of popula- 
tion differentiation for genetically neutral markers (fsx) 
and quantitative phenotypic traits (Qsx; Whitlock 2008) or 
on population genetics methods allowing to disentangle 
mean fitness estimates from the contribution of either gene 
flow and/or genetic drift (e.g., Whitlock and Gilbert 2012). 
Most of these models suggest that gene flow will tend to 
negate the selective effect and result in a lack of local adap- 
tation (or sometimes even in local maladaptation) in the 
pathogens, although many specific circumstances, related 
either to the genetics of the traits (dominance, epistasis) or 
to the spatial and temporal distributions of host and patho- 
gen populations can lead to apparent selection or to appar- 
ent neutrality on quantitative traits (Whitlock 2008; 
Blanquartet al. 2012). 

Unfortunately, there is a general lack of suitable experi- 
mental data to validate these theoretical predictions over a 
range of situations. Such experimental data would be best 
collected in situations (i) favoring metapopulation dynam- 
ics (geographically close but distinct sites with frequent 
local extinction) and (ii) with different dominant host 
genotypes between sites. This conjunction is well met in the 
potato late blight pathosystem close to the English Chan- 
nel. This pathosystem is indeed characterized by the follow- 
ing: 

A biotrophic pathogen 

The oomycete Phytophthora infestans - feeding on an 
annual host. The strict dependence of the pathogen upon 
living host tissue to persist in the environment and the 
annual rather than perennial growth habit of its hosts make 
periodic local extinction highly likely. This assumption is 
confirmed by genetic data, which showed that over 90% of 
the genotypes present in local populations in the Nether- 
lands disappear between one growing season and the next 
(Drenth et al. 1994); 
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A fragmented host distribution 

In western France, potato is commercially grown primarily 
in small coastal areas for the early table market. These spots 
are scattered along the northern coast of Brittany and of 
nearby Normandy. They are dedicated to vegetable produc- 
tion (typically potato, carrot, and leek in Normandy; 
potato, cauliflower, and artichoke in Brittany), with usually 
two different crops a year on each plot of land. This farther 
restricts the vegetation period of each of them and there- 
fore reinforces the chance for local extinction of biotrophic 
pathogens. Potato is a major crop in Jersey, one of the 
Channel Islands close to the coasts of Normandy and Brit- 
tany. 

Exposure to frequent and often strong winds 

The oceanic situation and mid-latitude of this part of Eur- 
ope goes along with a prevailing southwesterly wind 
regime, which may allow wind-borne spores to be dis- 
persed easUy over relatively large distances. 

Different dominant cultivars in Jersey and on the 
continent 

In the French potato production area considered, cultivar 
Europa was dominant but not exclusive. By contrast, 
potato production on Jersey has been dedicated almost 
exclusively since the early 1900s to cv. Jersey Royal, a selec- 
tion of the early cultivar International Kidney. Jersey Royal 
is grown year round in Jersey and is not produced any- 
where else. Jersey therefore represents a unique situation in 
Europe, because of its island location, the extensive mono- 
culture of a single potato cultivar for more than a century, 
and strict controls on movement of seed tubers of this cul- 
tivar into and out of Jersey (Deahl et al. 2009). Both culti- 
vars are susceptible to P. infestans, according to the 
European Potato Database (www.europotato.org). 

We used this unique situation to analyze experimentally 
the effect of gene flow and local selection on local adapta- 
tion in a host-pathogen system. Previous work showed that 
P. infestans populations from western France in 2004-2005 
involved the admixture of two genetically distinct groups 
of clonal lineages (Montarry et al. 2010a), both adapted to 
the prevalent cultivar in Northwestern Europe, Bintje 
(Montarry et al. 2008). Here, we wanted to determine 
whether gene flow exists between populations in Jersey and 
on the continent and would be sufficient to negate the 
effect of extensive, prolonged local selection by a single host 
cultivar, Jersey Royal. To this end, we tested the hypothesis 
of local adaptation to Jersey Royal in the P. infestans popu- 
lation present in Jersey, using both genotypic and pheno- 
typic data from isolates sampled in Jersey and in four 
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distinct French potato-cropping areas during a single epi- 
demic and a local adaptation cross-inoculation experiment. 

Materials and methods 

Collecting and maintaining Phytophthora infestans isolates 
During the summer of 2006, 232 P. infestans isolates were 
collected from 19 potato fields distributed across five sites 
in Jersey, Brittany (Ploudaniel, Paimpol and Saint-Malo), 
and Normandy (Val de Saire; Table 1). Fields within a site 
were separated by 1-10 km. At each site, samples were 
taken from the locally dominant cultivar: Jersey Royal in 
Jersey, Bintje in Ploudaniel (western Brittany), and Europa 
in the three other French sites (Paimpol, Saint-Malo and 
Val de Saire). Plants were haphazardly sampled across each 
field, and only one leaf per plant was sampled. Single lesion 
isolates were established and maintained as axenic cultures 
on pea agar medium as previously described (Montarry 
et al. 2010a). Isolates were kept at 15°C in the dark and 
transferred every 4-6 weeks to fresh agar media until phe- 
notypic characterization or DNA extraction. 

Mating type determination 

Phytophthora infestans is a heterothaHic species, with two 
mating types named Al and A2. The production of sexual 
offspring (oospores) occurs only when Al and A2 isolates 
come into contact. The mating type of each isolate was 
identified by pairing it on pea agar medium with known 
Al and A2 tester isolates, incubating in the dark at 15°C 
for 10-14 days, and looking under the microscope for the 
presence of oospores where the two colonies met. Isolates 
producing oospores with the Al tester were rated as A2 
mating type and vice versa. 

Genotyping isolates 

Isolates were grown in pea broth, previously autoclaved for 
20 min at 120°C and cooled to room temperature. After 
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10-15 days of incubation at 15°C, mycelium was washed 
three times in sterile water, dried on sterile paper, and 
lyophilized. DNA was extracted using the NucleoSpin® 
Tissue kit (Macherey-Nagel GmbH and Co, Dliren, Ger- 
many), according to the manufacturer's instructions. DNA 
concentration and purity were estimated for some isolates 
using a spectrofluorimeter Nanodrop ® ND-1000 (Fisher 
Scientific, IlUdrch, France). 

Fourteen polymorphic SSR loci: Pi02, Pi 04, Pil6, Pi33, 
Pi56, Pi63, Pi70, Pi89 (Lees et al. 2006), D13, Gil, 4B, 4G 
(Knapova and Gisi 2002), M19, and M34 (developed in 
our laboratory) were amplified using polymerase chain 
reaction (PGR). For Ml 9 and M34, respectively, the for- 
ward primers were 5'-cacgacgttgtaaaacgactagtggagacgttg 
TGCAGG-3' and 5'-cacgacgttgtaaaacgaccctttagcagctcatttc 
CG-3', and the reverse primers were S'-giggctgcggatacttctt 
T-3' and 5'-cgcactcgtcttgttgcac-3'. All PCRs were performed 
in a 12.5 /iL volume, as described by Montarry et al. (2010b). 
To detect simultaneously the alleles at several loci, the primers 
used were labeled with four fluorescent dyes: FAM (D13, 
Pi63, Pi02, Pi70, and M19), NED (Gil, 4G, Pi89, Pi33, and 
M34), VIC (Pi04 and Pil6), and PET (Pi56 and 4B). Nonla- 
beled primers were provided by Sigma-Genosys Ltd (Sigma- 
Aldrich, St. Louis, MO, USA), and labeled primers by AppHed 
Biosystems UK (Warrington, Cheshire, UK). The PCRs were 
performed in a MJ Research PTC-225 thermocycler under 
the following conditions: one cycle of 5 min at 94°C, fol- 
lowed by 30 cycles of 20 s at 94°C, 25 s at 58°C, and 30 s at 
72°C, and a final elongation cycle of 5 min at 72°C. 

Amplification products were pooled into four batches, 
based on expected allele sizes: batchl, Pi02, Pi89, and 4B; 
batch 2, Gil, Pi04, Pi56, Pi63, and Pi70; batch 3, D13, 4G, 
Pil6, and Pi33; batch 4, M19 and M34. For each batch and 
isolate, a 10-^L sample comprising 9.84 /iL deionized form- 
amide Hi-Di™ (Applied Biosystems), 0.06 /iL of 400 HD 
ROX™ size standards (Applied Biosystems), and 0.1 /iL of 
PGR multiplexed products were loaded into an ABI Prism 
3130x1 DNA Sequencer (Applied Biosystems) and run 
according to the manufacturer's instructions. DNA frag- 



Table 1. Sources of isolates of Phytophthora infestans collected in summer 2006. 

No. samples 

Location Potato cultivar Sampling date (No. fields) A1/A2 GPS central point 

Jersey Jersey Royal 15/05/2006 48 (5 fields) 48/0 49°13'02.25"N 

2°04'27.93"O 

Paimpol Europa 02/06/2006 46 (5 fields) 33/13 48''46'35.89"N 

3°04'34.95"O 

Saint-Malo Europa 22/06/2006 51 (2 fields) 51/0 48°37'0.45"N 

1°57'21;36"0 

Ploudaniel Bintje 06/07/2006 41 (3 fields) 40/1 48°30'06.10"N 

4°19,27.47"0 

Val de Saire Europa 12/07/2006 46 (4 fields) 46/0 49°39'21 .05"N 

1°16'38.44"0 
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ments were automatically sized with GeneMapper™ v4.0, 
and allele sizes were calibrated to the allele sizes of reference 
isolates (C1-C12) kindly provided by Drs A.K. Lees and 
D.E.L Cooke (James Hutton Institute, Dundee, UK). 

Quantitative pathogenicity 

Plants of the potato cultivars Bintje, Europa, and Jersey 
Royal were grown from certified seed tubers in a glasshouse 
regulated at 15-20°C (night/day temperatures) and 16 h of 
photoperiod, in 13-cm pots (one tuber per pot) filled with 
sand-peat-compost mixture. Plants were watered with a 
nutrient solution (Hakaphos; NPK 15/10/15) once a week. 
Leaflets were collected on 6- to 8-week-old plants. All 
experiments were performed from mid-September to mid- 
December 2007. 

To restore pathogenicity possibly lost during axenic cul- 
tures, each isolate was first multiplied separately on 
detached leaflets of the susceptible cultivar Bintje. Each 
leaflet was placed abaxial side up and infected by deposing 
a drop of a suspension of P. infestans sporangia, prepared 
by flooding a 3-week-old culture with 5 mL sterile distilled 
water and gently scrapping the colony surface with a small 
glass rod. After 7 days of incubation in humid chambers 
under controlled conditions (15°C/18°C night/day temper- 
atures, 16-h daylight), the new sporangia produced on the 
infected leaflets were collected in sterile water, and the con- 
centrations of the resulting suspensions were adjusted to 
5 X 10* sporangia mL^'. Suspensions were kept at 4°C for 
approximately 2-3 h to promote zoospore release, just 
before being used for pathogenicity experiments (Montarry 
et al. 2008). 

Each isolate was tested for lesion size and spore produc- 
tion on aU three host cultivars (Bintje, Europa, and Jersey 
Royal). Six leaflets per cultivar x isolate combination were 
placed abaxial face up on the lids of inverted Petri dishes 
containing 10 g L^' water agar (two leaflets per dish) and 
infected by depositing a 20 /(L drop of the calibrated spo- 
rangial suspension at the leaflet center. Infected leaflets 
were incubated for 6 days as described above. Lesion diam- 
eters were measured in two perpendicular directions, and 
lesion size was computed assuming an elliptic shape. Then, 
each leaflet was washed in 10 mL saline buffer (Isoton II), 
and sporangia were counted with a Beckman Coulter Z2 
cell counter (Beckman, Villepinte, France) to determine 
sporangial production per lesion. 

Data analyses 

Genetic data 

Among the fourteen microsatellite loci, only twelve were 
analyzed because over 15% of the isolates showed multiple 
banding at loci M19 (28 isolates showing three bands and 
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five showing four bands) and Pi63 (all isolates showing 
three bands), casting doubt as to the actual ploidy level at 
these loci. Moreover, to have a correct assignment to mul- 
tilocus genotypes, all isolates with amplification in which 
amplification failed for at least one locus were excluded 
from the final data set, which therefore included 179 P. 
infestans isolates without any missing data. 

Repeated multilocus genotypes (MLGs), that is, isolates 
sharing the same alleles at all 12 loci, and the number of 
copies of each MLG (n) were identified using GenClone 
2.0 (Arnaud-Haond and Belkhir 2007). Genotypic diver- 
sity was estimated as the genotypic richness R = (num- 
ber of different MLGs - l)/(total number of isolates — 1). 
Because an MLG may result from distinct sexual repro- 
duction events or from clonal reproduction, we estimated 
the probability that the n occurrences of a given MLG 
results from sexual reproduction, given the observed allele 
frequencies for each population (Psex n re-encounter). 
Because the inclusion of clonal multicopies can strongly 
distort the _F-statistics, we considered only one copy of 
each multilocus genotype for aU subsequent genetic analy- 
ses {i.e., a clone-corrected data set). An exception is the 
computation of the genetic diversity indices, for which 
statistics were performed both with and without clonal 
copies. 

The evolutionary relationships among P. infestans MLGs 
were inferred from a minimum spanning network, calcu- 
lated using MiNSPNET (Excoffier and Smouse 1994) from the 
matrix of genetic distances based on the number of differ- 
ent alleles (DAS), and represented graphically via the neato 
layout engine using graphviz (http://graphviz.org). 

Linkage disequilibrium between loci was calculated with 
GENEPOP (with default Markov chain parameters) to count 
the number of locus pairs showing significant linkage dis- 
equilibrium across all populations. A Bonferroni correction 
(adjusted a = 0.0009; 55 comparisons) was applied to take 
into account multiple testing. 

Genetic diversity in each population (and each cluster) 
was estimated as both the allelic richness (Ar) and gene 
diversity. Allelic richness was estimated using the rarefac- 
tion method implemented in populations 1.2.30 (LangeUa 
1999), which estimated the mean number of alleles per 
locus for a reduced sample size. An unbiased estimate of 
gene diversity (Hg according to Nei, 1978) and deviation 
from random mating (Fis) were computed using genetix 
4.05.2 (Belkhir et al. 1996-2004). The statistical signifi- 
cance of _Fis values for each population (and each cluster) 
was tested using the resampling method (1000 bootstraps) 
implemented in genetix. 

The differentiation coefficients between each pair of 
P. infestans populations (-Fst) were estimated according 
to Weir and Cockerham (1984) using genepop 4.0.7 (Ray- 
mond and Rousset 1995), and their statistical significance 
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was tested with the exact G-test (with default Markov 
chain parameters) implemented in genepop. Following the 
recommendations of Rousset (1997), the hypothesis of 
isolation by distance was tested by calculating the correla- 
tion between the matrices representing pairwise Slatkin 
genetic distances Fsrlil—FsY, Slatkin 1995) and the natu- 
ral logarithm of geographic distance for each pair of 
P. infestans populations. The statistical significance of the 
correlation was assessed with a Mantel test (10 000 per- 
mutations) using XLSTAT 2012 (Addinsoft SARL, Paris, 
France). 

Without taking into account their geographic origin, 
P. infestans isolates were clustered on the basis of their 
genetic relatedness using the Bayesian clustering approach 
implemented in structure 2.3.3 (Pritchard et al. 2000; 
Falush et al. 2003). Simulations were performed using the 
admixture model. We estimated the number K of genetic 
clusters (here between K = 1 and K = 15) to which the iso- 
lates should be assigned. For all simulations, we did not 
force the model with predefined allele frequencies for 
source clusters. Five independent runs were conducted to 
assess the consistency of the results across runs, and all runs 
were based on 500 000 iterations after a burn-in period of 
100 000 iterations. We then identified the number of 
genetically homogeneous clusters as described by Evanno 
et al. (2005). 

Phenotypic data 

Normality and homogeneity of variances were checked 
with the Shapiro-Wilk and the Leven tests, respectively, 
using the statistical software R, version 2.15.0 (The R Foun- 
dation for Statistical Computing, 2012). The effects of 
potato cultivar, P. infestans population, and the corre- 
sponding interaction on lesion size and spore production 
were tested through an anova model using the statistical 
software R. A second anova model was used to test the 
P. infestans population effects separately for each potato 
cultivar. Indeed, when environments are defined as differ- 
ent host resistance levels, the 'home versus away' criterion 
is poorly informative regarding adaptation to the host, 
because it is biased by an important host effect, and the rel- 
evant test for local adaptation becomes the 'local versus 
foreign' criterion (Thrall et al. 2002; Kawecki and Ebert 
2004; Montarry et al. 2008). When significant effects were 
detected, mean values were compared with Tukey's HSD 
tests (a = 0.05). 

Results 

Mating type 

Most isolates belonged to the Al mating type: only 14 iso- 
lates of the 232 isolates were A2 and were found in one field 
sample at Paimpol (13 isolates) and one isolate in Ploudan- 
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iel (Table 1). Among those 14 A2 isolates, 12 isolates, all 
from Paimpol, were successfully genotyped. 

Multilocus genotype analysis 

Forty-four alleles could be amplified across the 12 microsat- 
ellite loci, with two to seven alleles per locus. GenClone 
identified 23 unique multilocus genotypes (MLGs) among 
the 179 P. infestans isolates, corresponding to a genetic 
richness of _R = 0.128. Fifteen MLGs were represented by a 
single isolate (Table 2). Three of the eight repeated MLGs 
were found in at least two populations: (i) MLG9 included 
four isolates from Paimpol and nine isolates from Ploudan- 
iel; (ii) MLG16 included 15 isolates from Jersey, seven iso- 
lates from Paimpol, 10 isolates from Saint-Malo and four 
isolates from Val de Saire; and (iii) MLG18 included 19 iso- 
lates from Jersey, 12 isolates from Paimpol, 20 isolates from 
Saint-Malo, 11 isolates from Ploudaniel, and 31 isolates 
from Val de Saire. The very low P^^x n re-encounter values 
indicated that the overrepresentation of these repeated 
MLGs probably results from clonal amplification (Table 2). 

Two MLGs (MLGl and MLG2) included only isolates 
belonging to the A2 mating type (Table 2). Comparing 
their fingerprints with those described by Cooke et al. 
(2012) showed that they differed markedly from the domi- 

Table 2. Characteristics of the 23 multilocus genotypes (MLGs) which 
were discriminated using 12 microsatellite loci on the 179 Phytophthora 
infestans isolates. Cluster assignment, number of copies (n), mating 
type, and Pjex n re-encounter are indicated for each MLG. 





Cluster 








MLG 


assignment 


n 


Mating type 


re-encounter 


MLG 18 




93 


Al 


2.67E-256 


MLGl 6 




36 


Al 


6.50E-110 


MLGl 2 




6 


Al 


2.31E-1B 


MLG21 




2 


Al 


4.18E-0B 


MLG22 




2 


Al 


1.53E-06 


MLG5 






Al 




MLG6 






Al 




MLG11 






Al 




MLG 13 






Al 




MLG 14 






Al 




MLG 15 






Al 




MLG 17 






Al 




MLG 19 






Al 




MLG20 






Al 




MLG23 






Al 




MLG9 


2 


13 


Al 


2.10E-84 


MLG2 


2 


9 


A2 


4.24E-48 


MLGl 


2 


3 


A2 


2.25E-13 


MLG3 


2 




Al 




MLG4 


2 




Al 




MLG7 


2 




Al 




MLGS 


2 




Al 




MLG10 


2 




Al 
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nant MLGs in Great Britain at the time, that is, 13_A2, and 
were rather variants of MLG 3_A2, a rapidly decHning hne- 
age in the British Isles between 2005 and 2007. A similar 
observation was made regarding the main Al MLGs, which 
were variants of either 1_A1 (MLG9) or 2_A1 (MLG12, 
-16, and -18) rather than of either 6_A1 or 8_A1, which 
were the major Al lineages in Britain at the time. MLG18, 
which is the most abundant within our samples, is also 
similar to MLG JE-1.5 reported in Jersey during the same 
period (Deahl et al 2009). 



Table 4. Pairwise F^j distances between Phytophthora infestans popu- 
lations. 





Jersey 


Val deSaire 


Saint-Malo Paimpol 


Val deSaire 


-0.012 






Saint-IVlalo 


-0.023 


-0.008 




Paimpol 


0.009 


0.062* 


0.035 


Ploudaniel 


0.152* 


* 0.194*** 


0.182*** 0.023 


Significance 


of FsT is 


indicated by stars 


(*<0.05, **<0.01 and 



***<0.001). 



Genetic structure of Phytophthora infestans populations 

Four to nine MLGs were found in each of the five popula- 
tions, the population from Ploudaniel showing the highest 
genotypic diversity [R = 0.30; Table 3). Whatever the data 
set used (with or without clonal copies), the mean allelic rich- 
ness (Ar) per population ranged from 1.88 to 2.41 and gene 
diversity (He) from 0.383 to 0.514 (Table 3). The highest Ar 
and He were found in populations from Paimpol and Plou- 
daniel. Fis estimates for each population were in the majority 
of cases significantly negative, ranging from —0.769 to 
—0.151 and from —0.645 to 0.018 for data sets with and with- 
out clonal copies, respectively (Table 3). Moreover, 12 locus 
pairs among the 55 possible comparisons showed significant 
linkage disequiHbrium (P < 0.05) across aU populations if no 
corrections for multiple tests across all populations were 
made, although none showed significant linkage disequilib- 
rium after Bonferroni's correction. 

The genetic differentiation between populations was sig- 
nificant in four of the ten pairwise Pst values (Table 4). 
The population from Ploudaniel was significantly highly 



differentiated {P < 0.01; Table 4) from those from Val de 
Saire (Pst = 0.194), Saint-Malo (Psr = 0.182), and Jersey 
{FsT = 0.152). To a lesser extent, populations from Paim- 
pol and Val de Saire were also genetically differentiated 
{FsT = 0.062). Consequently, a significant correlation was 
observed between geographic distance and P. infestans 
population differentiation (Mantel test, « = 9; Pearson's 
coefficient of correlation r = 0.87; P = 0.0009; Fig. 1). The 
isolation by distance pattern was, however, nonsignificant 
when the population from Ploudaniel was excluded 
(r = 0.65, P = 0.164), consistent with the low differentia- 
tion among the easternmost populations sampled (Table 4, 
Fig. 3). 

Clustering analysis and geographic repartition of genetic 
clusters 

A Bayesian clustering analysis was performed using struc- 
ture on the clone-corrected data set {i.e., 23 P. infestans 
MLGs). That clustering analysis clearly indicated that the 
posterior distribution of the allele frequencies among clus- 



Table 3. Genotypic and genetic diversity indices for each Phytophthora infestans population (A) and for eacin cluster (B). N, number of isolates; MLG, 
number of distinct multilocus genotypes; R, genotypic richness ([\/1LG-1)/(N-1); He, unbiased estimate of gene diversity (Nei, 1978); Ar, allelic richness 
corrected for sample size (n); f,; significantly different to zero are indicated in bold. The genetic statistics (Hj, Ar, and F,;) were performed both with 
and without clonal copies. 

Genetic diversity 





Genotypic diversity 




With clonal 


copies 






Without clonal copies 




A 


W 


MLG 


R 


He 


Ar (n = 


= 28) 


fis 


He 


Ar (n = 4) 


Fis 


Jersey 


36 


4 


0.09 


0.384 


1.88 




-0.765 


0.470 


1.92 


-0.645 


Paimpol 


36 


6 


0.14 


0.462 


2.41 




-0.369 


0.514 


2.33 


-0.212 


Saint-Malo 


38 


5 


0.11 


0.383 


1.91 




-0.757 


0.456 


1.91 


-0.553 


Ploudaniel 


28 


9 


0.30 


0.433 


2.33 




-0.151 


0.443 


2.19 


0.018 


Val de Saire 


41 


7 


0.15 


0.378 


1.95 




-0.769 


0.449 


1.93 


-0.486 










Genetic diversity 














Genotypic diversity 




With clonal copies 






Without clonal copies 




B 


W 


MLG 


R 


He 


Ar (n = 


30) 


Fis 


He 


Ar (n = 4) 


Fis 


Cluster 1 


149 


15 


0.09 


0.381 


2.06 




-0.751 


0.459 


2.25 


-0.451 


Cluster 2 


30 


8 


0.24 


0.333 


2.08 




-0.422 


0.344 


2.08 


-0.264 
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Figure 1 Isolation by distance pattern between genetic differentiation, 
measured as fsiAl - Fst). and geographic distance (natural logarithm 
of the distance in km) for pairwise Phytophthora infestans populations. 



ters was best explained with a grouping into K = 2 genetic 
clusters. Applying an assignment threshold of 0.90, 91% of 
the MLGs (21 of 23) could be assigned to one or the other 
genetic cluster, confirming that the two clusters were well 
differentiated (Fsi = 0.317). The evolutionary relationships 
among the P. infestans MLGs from both genetic clusters 
were showed by a minimum spanning network (Fig. 2). 

Cluster 2 had a higher genotypic diversity (R = 0.24 with 
eight MLGs for 30 isolates) than cluster 1 (R = 0.09 with 
15 MLGs for 149 isolates; Table 3). Although cluster 1 
included only isolates of the Al mating type, whereas clus- 
ter 2 included both mating types (Al and A2), and Ar 
were higher for cluster 1 than for cluster 2, especially in 
analyses performed without clonal copies (Table 3). Both 
clusters showed significantly negative Fis (Table 3). While 
cluster 1 included isolates from all populations, only iso- 
lates from Ploudaniel and Paimpol were present in cluster 
2 (Fig. 3). 

Pathogenicity of Phytophthora infestans populations on 
different cultivars 

All tested explanatory factors [i.e., P. infestans populations, 
potato cultivars, and their interaction) had a significant 
effect on both lesion size and spore production (Table 5). 
Overall, Jersey Royal cv. tended to be less susceptible than 
the other two cultivars (Fig. 4), and the population from 
Jersey tended to be less pathogenic than the French ones, 
irrespective of the cultivar tested (Fig. 4), explaining the 
significant population effect in the anova (Table 5) and 
suggesting general maladaptation in this population. The 
second anova model showed significant P. infestans popu- 
lation effects for each potato cultivar (Bintje, Europa, and 
Jersey Royal) and each pathogenicity trait (lesion size and 
spore production), except for lesion size on Jersey Royal 
where the population effect was only marginally significant 
(f3_i57 = 2.54, P = 0.058). 



Discussion 

Earlier work with geographically distant populations has 
shown that local adaptation can emerge in P. infestans 
(Andrivon et al. 2007), but is negated when populations 
are strongly connected and/or submitted to largely uniform 
selection pressures by a single host cultivar (Montarry et al. 
2008). Here, we expected local adaptation to be detectable 
in Jersey, because of its insularity and of prolonged selec- 
tion by a host cultivar specific to this island. Furthermore, 
the mUd climate and year-round presence of the host were 
expected to limit the extinction probability of locally 
adapted genotypes and therefore to reinforce such a local 
adaptation pattern. 

However, our data revealed no pattern of local adapta- 
tion to cv. Jersey Royal in Jersey, and extensive gene flow 
between P. infestans populations in Jersey and in the 
neighboring parts of France. The low genetic differentia- 
tion between the population from Jersey and its three 
nearest neighbors from France (Val de Saire, Saint-Malo 
and Paimpol) supports the hypothesis of intensive gene 
flow at this scale, despite the apparent phenotypic differ- 
ences in pathogenicity. Together, these observations 
strongly suggest that gene flow between France and Jersey 
is preventing the evolution of local adaptation patterns, 
despite persistent and widespread selection in Jersey over 
more than a century by a single host cultivar specific to 
this island. 

The geographic differentiation between local populations 
gives clues as to the effective distance upon which gene flow 
is active in this pathosystem. This distance is in the order 
of 100-120 km, according to the F^i calculations between 
pairs of populations. It explains the significant isolation by 
distance pattern when the Ploudaniel population is consid- 
ered, but the lack of significance of this pattern is explained 
when this population, further west, is excluded from the 
data set. The dispersal range calculated here markedly 
exceeds the usual assumptions, of the order of magnitude 
of a few kms, regarding the effective dispersal kernels of the 
short-lived, fragile sporangia of P. infestans (see Skelsey 
et al. 2010). This may be due to the very conducive climatic 
conditions favoring dispersal of the pathogen in western 
France, with strong winds, widespread potato crops, and 
frequent episodes of high humidity, which all favor infec- 
tion by disseminated spores. 

A further factor favoring effective dispersal and general 
adaptation is the genetic proximity between hosts. Search- 
ing the European Potato Database (www.europotato.org) 
reveals that cv. Europa has cv. Bintje among its progenitors, 
through CVS. Climax and Edzina. It is therefore understand- 
able that adaptation to cv. Bintje, previously shown in 
French populations (Montarry et al. 2008), might also con- 
fer adaptation to some of its offspring. This then favors the 
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Figure 2 Minimum-spanning network sliow- 
ing tine relationships among the Phytophthora 
infestans multilocus genotypes (MLGs 1 to 23) 
of both clusters (1 in white and 2 in gray). 
Branch sizes are proportional to genetic dis- 
tance {/. e., the number of different alleles) and 
circle areas to the numbers of isolates. 

persistence of a large population with the same core adap- 
tive features and thus limits the possibility for emergence 
of local adaptation patterns and metapopulation dynamics 
with divergent selection histories. 

Our finding that P. infestans populations cluster into two 
distinct genetic groups is consistent with previous data from 



French isolates collected in 2004 and 2005 (Montarry et al. 
2010b) revealing that the admixture pattern was stable over 
at least three successive years (from 2004 to 2006). The 
maintenance of the differentiation between those clusters 
could be explained by the fact that sexual reproduction 
seems to play no role in the adaptive patterns recorded here. 
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Figure 3 Geographic position of each Phytophthora infestans popula- 
tion used in this study and mean assignment percentage for each popu- 
lation to each genetic clusters. The populations linked by black lines are 
not genetically differentiated (P > 0.05) based on an Fsj test, and those 
not linked by black lines are genetically differentiated (P < 0.05). 

The population structures observed remained highly clonal, 
although admixed, and no sign of recombination was 
observed despite the co-occurrence of both mating types in 
two of the populations surveyed (Paimpol and Ploudaniel). 
Indeed, the P^ex values and the negative f is values, indicat- 
ing an excess of heterozygotes that could result from asexual 
reproduction (Halkett et al. 2005; Goyeau et al. 2007), 
bring strong evidence for clonal reproduction in these pop- 
ulations. This general structure, consistent with earlier 
observations (Montarry et al. 2010b), probably further 
strengthens the persistence of general adaptation to Bintje 
rather than the emergence of local adaptation. 

These observations have major consequences for resis- 
tance deployment strategies and the management of late 
blight epidemics through resistant cultivars. Such deploy- 
ment strategies, to be effective, should ideally exploit local 
maladaptation generated when the hosts are diverse and 
the pathogen moves rapidly (Kawecki and Ebert 2004). 
They can also rely on a strong metapopulation dynamics 
resulting from limited gene flow and frequent local extinc- 




Bintje Europa Jersey Royal 

Potato cultivar 



30 




Bintje Europa Jersey Royal 

Potato cultivar 



Figure 4 Mean lesion size and spore production of each Phytophthora 
infestans population, from Jersey, Paimpol, Saint-Malo, and Val de 
Saire, on the different potato cultivars, Bintje, Europa, and Jersey Royal. 
Letters represent the homogenous groups identified with the Tukey's 
HSD test at the 5% threshold. 

tion. The strictly biotrophic nature of P. infestans and its 
host range limited to a few annual hosts should favor such 
local extinction. However, the widespread presence of 
potato crops throughout the territory and the extensive 
gene flow shown here imply that a metapopulation dynam- 
ics is quite unlikely to take place. This is consistent with 
population structures of P. infestans in western Europe 
being dominated by a few, widely distributed clonal lin- 
eages, such as US-1 until the early 1980s (Spielman et al. 
1991) or 13_A2 recently (Cooke et al. 2012). 



Table 5. Anova of lesion size and of spore production. Sources of variation are Phytophthora infestans population, potato cultivar, and the corre- 
sponding two-way interaction between these variables. 

Lesion size Spore production 

Mean Mean 
Source of variation df square f value P>F df square Fvalue P>F 



Phytophthora infestans population 3 

Potato cultivar 2 

Population * cultivar 6 

Error 478 



1212.6 51.89 <0.0001*** 3 
318.7 20.46 <0.0001*** 2 
308.2 6.59 <0.0001*** 6 

3723.7 472 



5.0E + 11 28.17 <0.0001*** 

2.9E + 11 24.82 <0.0001*** 

1.9E + 11 5.51 <0.0001*** 
2.8E + 12 



Statistically significant effects are indicated by asterisks (*P < 0.05; **P < 0.01; ***, P < 0.001); df, degrees of freedom. 
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Therefore, it will probably be difficult to think of effec- 
tive resistance deployment strategies over large scales, such 
as crop mosaics or refuge areas, because the strong gene 
flow and connection of host populations will make it hard 
to generate sufficient differentiation in local parasite pop- 
ulations for a large-scale spatial deployment strategy of 
resistant hosts to be effective (Skelsey et al. 2010). The use 
of highly differentiated host genetic backgrounds, for 
instance the introduction of different resistant R gene pyr- 
amids through ad-hoc cassettes in a cisgenesis transforma- 
tion process (Haverkort et al. 2008), might help, although 
the durability of gene pyramids has proven low in the past 
due to the large effective population sizes of the pathogen 
(e.g., Andrivon 1994). A better option might thus be to 
promote within-crop resistance diversification through 
cultivar or species mixtures (Garrett and Mundt 2000; 
Andrivon et al. 2003; Philips et al. 2005; Bouws and Fin- 
ckh 2008), which minimizes genotype unit area and has 
proven effective when disease pressure can be maintained 
at a moderate level (Pilet et al. 2006), rather than diversi- 
fication strategies targeted at the landscape or regional 
spatial scales. 
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